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We present the properties of a cosmic superstring network in the scenario of flux compactification. 
An infinite family of strings, the (p,q)-strings, are allowed to exist. The flux compactification 
leads to a string tension that is periodic in p. Monopoles, appearing here as beads on a string, 
are formed in certain interactions in such networks. This allows bare strings to become cosmic 
necklaces. We study network evolution in this scenario, outlining what conditions are necessary 
to reach a cosmologically viable scaling solution. We also analyze the physics of the beads on a 
cosmic necklace, and present general conditions for which they will be cosmologically safe, leaving 
the network's scaling undisturbed. In particular, we find that a large average loop size is sufficient 
for the beads to be cosmologically safe. Finally, we argue that loop formation will promote a scaling 
solution for the interbead distance in some situations. 
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I. INTRODUCTION 

Braneworld models of inflation can produce cosmic 
strings [H, [2, [2] • This discovery has reignited interest in 
cosmic strings as cosmological observables (see, e.g., 0,[1] 
for reviews). Since cosmic strings can also be formed in 
ordinary gauge-theory inflationary models, string theo- 
rists are eager to find ways in which the cosmic "su- 
perstrings" from brane inflation are different from the 
Abelian Higgs cosmic strings typically found in gauge- 
theory inflation. 

Among these distinctive phenomena, the best known 
is that higher-dimensional effects can cause the super- 
string intercommutation probability to be considerably 
less than one [2, @|; this reduced probability does not 
significantly alter network dynamics in most cases Q. 
However, an even more exotic new property of cosmic 
superstrings is that they can come in more than one fun- 
damental species. The braneworld permits the forma- 
tions of both fundamental cosmic strings, or F-strings, 
and Dirichlet brane cosmic strings, or D-strings. This 
leads to a new type of string interactions: binding of 
multiple strings. The network-level scaling behavior of 
these multi-tension, binding strings was first studied in 
Ref. Q. In that work, however, they assumed a simple 
string tension formula that is only valid for strings living 
in a flat 10 dimensional Minkowski space, which is not a 
realistic scenario. An important goal of this paper is to 
establish conditions under which cosmic superstring net- 
works can reach a scaling solution when they are realized 
in a more realistic - i.e., a phenomenologically viable - 
background geometry. 

Constructing more-realistic string theory models of the 
early universe has led to a great deal of work on how 
the complexities of higher-dimensional geometry are tied 
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to physics. In particular, braneworld models of inflation 
have been realized in a geo metry known as the Klebanov- 
Strassler (KS) throat flOf. Inspired by these develop- 
ments, Firouzjahi et al. [ll| showed that the properties 
of cosmic superstrings can depend on the structure of this 
extra-dimensional geometry. 

It has also been discovered that fundamental strings 
in a KS throat are dual to confining flux tubes. Since 
confining flux tubes are important in many contexts, dis- 
covery of this relationship has occasioned yet more inves- 
tigation of these strings. In Ref. [12], it was found that 
fundamental strings can join on a D3-brane wrapping the 
non-vanishing S 3 cycle that lives at the bottom of a KS 
throat. Since this is a three-dimensional object being 
wrapped on a three-dimensional cycle, the wrapped D3- 
brane looks like a point charge, or a monopole, in 4D. 
In the dual gauge theory, this monopole is equivalent to 
a baryon. Thus, cosmic superstrings networks can also 
contain monopoles. Since these monopoles are dual to 
baryons in the gauge theory, one might call these stringy 
monopoles "baryons," too. However, since this would 
likely cause confusion, we will instead call these objects 
"beads." These "beads" are stuck on cosmic strings and 
will act like beads on the cosmic "string" . Cosmic strings 
from string theory are thus actually cosmic "necklaces." 
Such string-bead "necklace" objects were first studied in 
the context of gauge theory defect models in Ref. [12]. 
Also, a different approach to cosmic necklaces in string 
theory can be found in Refs. [l4l [TEl]. 

As it turns out, Fundamental strings living in the 
warped KS throat geometry have a tension spectrum that 
is sinusoidal in its dependence on winding number, or flux 
quantum number. The period is set by the amount of RR 
F3 flux in the extra-dimensions. This is an integer de- 
noted M and the fact that the tension is periodic is an 
important new property that we must take into account. 
In this paper, we will treat M as a free parameter. Wc 
will show that requiring safe cosmological behavior al- 
ready places lower limits on what values M can have. 
To summarize: a collection of M F-strings can join on a 
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bead. They are Zm charged under the usual NS-NS 2- 
form J5 2 (which in 4D is Hodge dual to an axion). These 
facts, that cosmic F-strings are Zm charged and can join 
on beads, were first discussed by Ref. [la ]. 

Dirichlet strings, or D-strings, on the other hand, are 
believed to be BPS in the KS throat [izj]. The F-strings 
and D-strings together form a (p, q) string network. This 
network is so named because the mutual interactions be- 
tween F and D strings are binding interactions. Thus, a 
series of interactions between p F-strings and q D-strings 
can result in the formation of a bound state contain- 
ing the net charge for all the bound strings. This string 
bound state is called a (p, q) string. The tensions of the 
various (p, q) states in a warped throat were calculated 
in Ref. 0] to be 
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where m s — y is the warped string scale, g s is the 
string coupling constant and b is a numerical factor ap- 
proximately equal to 0.93. This formula was also derived 
using different methods in Refs. (l9l. [20j. 

As we noted previously, multi-tension networks of 
binding (p, q) strings have been studied before 6, 9]. One 
goal of our present work is to generalize the work of Ref. 
9] using the more realistic tension spectrum (Eqn. [1]) and 
taking into account periodicity in the quantum number, 
p. We also hope to understand better the "beads" that 
can be formed dynamically on these strings. For future 
reference, the mass of the beads is given by [18| 
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This paper is roughly divided into two different parts. 
The first is concerned with the scaling of the network 
while neglecting the beads. The second part discusses the 
physics of the beads themselves. We study the scaling of 
the network neglecting the beads using a modified version 
of the VOS model [2l[ . As this study will closely follow 
that of Ref. @ , we will not provide a detailed review of 
that model. 

The main result for the first part of our paper is that 
the network does have a scaling solution when M > 10. 
The more complex background geometry also leads to 
new network features. For instance, the structure of the 
tension spectrum is much richer than when using the 
flat space formula for the tension, and the overall cos- 
mic string density is slightly higher. Interestingly, we 
also find that the scaling density of cosmic strings de- 
pends on the network's initial conditions, with there be- 
ing, in effect, two possible scaling solutions. If the large 
quantum number string states - those with p ~ M - are 
excited initially, then the network will scale with com- 
parable number density of (1,0) and (M — 1,0) strings, 
effectively doubling the scaling network energy density as 
compared with the flat space case. We remind the reader 



that (1,0) and (M — 1,0) strings have the same tension. 
On the other hand, if the post-inflationary network cre- 
ation excites exclusively low quantum number strings - 
those with p < M/2 - then the network will reach scaling 
at a density only slightly higher than the flat space case. 
Though the former scenario is intriguing, we consider the 
latter to be more probable. 

All of the results summarized here were derived under 
the assumption that the beads' contribution to the net- 
work's energy density is negligible. We examine this issue 
in detail in section HVl To ease discussion of this issue, 
we will follow Berezinsky and Vilenkin [22| by quantify- 
ing the beads' contribution to the energy density with a 
dimensionless parameter, r: 
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where d is the average interbead distance along the string 
and /j, is the tension of the string. The energy per unit 
of length of a bead-carrying string is, then, (r + l)/x. 
For small r, the beads contribute a negligible amount 
to the energy density and can be neglected, while for 
large r they are cosmologically disastrous: the network 
becomes a gas of monopoles connected by relatively light, 
irrelevant strings. 

We first analyze the evolution of r in the scaling 
regime, neglecting loop formation, annihilation of beads, 
and taking the assumption that the network reaches scal- 
ing with a small value of r after an early, transient regime. 
In this case, we find that r grows unless the average size of 
loops (which controls the rate of gravitational wave emis- 
sion and thus string contraction) is bigger than about 
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where t is the FRW coordinate of time. 

Loop formation should, in general, make r increase as 
well, even if the average loop size respects the previous 
bounds. But taking account of loops will not always force 
r to grow indefinitely Instead, wc believe r should reach 
a scaling value, with the interbead distance of the same 
order as the average loop size. We argue for this view 
more extensively in section IIVB1 In this approach, the 
approximate scaling value for r should be given by 



10 2 M 3/: 



V2glb H 
37r 3 / 2 m, 



(5) 



where H ~ 1/t is the Hubble constant. Strangely, un- 
der these assumptions we find that larger values of the 
flux, M, lead to larger values of r and, thus, potentially 
more observable features. This is a counter-intuitive re- 
sult, since in the large M limit we recover the flat space 
tension formula, which does not share any of these new 
features. Dynamical bead formation is also heavily sup- 
pressed for large values of M. Thus, we should regard this 
result as suggestive but incomplete. The proportionality 
between r and M reflects the fact that monopole mass 
grows with M more quickly than does the string tension. 
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An interim regime, where r ~ 1, would be very inter- 
esting to study. In this case, the effect of beads would 
cosmologically be safe, yet observable: essentially dou- 
bling the effective tension of every string in the network. 
However, even in the most optimistic case, the fact that 
r oc H means that any period of bead-string energy bal- 
ance will be short lived; as H shrinks with cosmological 
expansion, r will also shrink. It is thus unlikely that there 
will be a great number of necklace-specific effects that re- 
main important long enough to be easily observable. 

Finally we comment on the physics of the beads in the 
transient regime, when the network evolves from its ini- 
tial configuration to a scaling solution. This is a much 
murkier issue. Only a complete simulation can really ad- 
dress this epoch with any certainty. Yet it is possible 
to get an idea for what may happen during this regime 
by employing the same logic we use in the scaling era. 
We expect r to increase significantly during this regime 

during which we expect a great number of string inter- 
actions to occur - so that we will need a large value of 
AI in order to prevent r from growing beyond 1. If we 
take quantum mechanical tunneling suppression of bead 
formation as the main mean for limiting the growth of 
r, we find that M should be ~ 10 3 or more to control 
the growth of r during the transient regime. But, again, 
this result should be taken as very conjectural, given the 
complexity of network physics outside of the scaling era. 

Each of these issues will be better fleshed out below. 
In section UH we summarize the main properties of the 
cosmic string network. In section InTI we discuss the net- 
work's scaling solution, obtained through numerical in- 
tegration of a modified MTVOS network model. Section 
|IV] treats the physics of the beads in more detail. 



II. INTERACTIONS 

We expect no substantial change from the calculation 
of Ref. 6] for local binding events. Differences appear 
only when the interaction forms beads. 



A. Probability of Interaction 

As in Ref. 0], when two strings of different types in- 
teract they form a Y-junction. Defining the quantum 
numbers p and q to be positive for a string that is di- 
rected outward from the junction, charge conservation 
implies that: 



Y,p = nM > 
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where n, the bead number, is valued in Z (negative n 
implies an antibead). Thus, any interaction that results 
in p > M implies that a bead - antibead pair must be 
formed to conserve charge (cf. Fig. [I} . In a fully general 
treatment, we should note, second order flux effects will 
come into play, complicating this story (and lifting the 



tension of the M-string above zero) , though we expect no 
qualitative changes to be wrought by these considerations 
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FIG. 1: This shows the intercommutation of a (M-1,0) string 
and a (1,0) string. The strings have the same tension but dif- 
ferent charge. They form a Y-junction with an M-string with 
zero tension, which in some sense invisibly connects the bead 
to the anti-bead. The beads are the finite energy objects that 
mediate the change of topological charge from one value to 
the other. One can have beads and antibeads (conventionally 
chosen red/empty and blue/filled on the picture). We also 
show the formation of a closed loop via this process. 

We can calculate the probability of bead formation 
by treating the formation of a bead-antibead pair as a 
tunneling event. Consider a simple triangular potential 
V(x), with V(0) = 2M b and V(L) = for L = 2M b /fi. 
Using the WKB approximation, the probability of inter- 
action is P ~ e~ 7rA , where A is the area under the curve. 
The factor of n is included to agree with a more precise 
calculation using an instanton method. We then find that 



tunneling ~ e 
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Interestingly, for the case of Fig. [TJ the probability is 
close to one for an interesting range of M values 
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We see that A. — ► oo as M — ► oo, as we should expect, 
since the bead mass grows faster than the string tension 
as a function of M. However, since the factor that multi- 
plies M 2 may be of order 10 -3 (if we assume g s ~ 10 _1 ), 
we should have a non-suppressed production of beads for 
M < 10- Note that for smaller g s this minimum value 
will increase. 

In Ref. [Hj], they argue that M needs to be greater 
than 12 in order for any anti-D3-brane to be classically 
stable in the KS throat. Furthermore, we need M to be 
greater than ~ l/Ss in order for the size of the extra 
dimensions to be big enough to trust the supergravity 
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approximation. These considerations imply that bead 
production is always suppressed, albeit only slightly if 
M is not much above these bounds. 

When two strings of different tensions meet and form 
beads (see Fig. [2|, there are two different probabilities 
for tunneling since we have two different incoming string 
tension from which to choose. When this occurs in our 
calculations, we calculate the probability using the larger 
tension, since this will be the most probable outcome. So 
we use the same formula as before, Eqn. [JJ with /i being 
the tension of the heavier of the two colliding strings. 

Finally, we have the possibility of multiple beads get- 
ting stacked on each other through network dynamics. 
Though we expect this to be rare if beads are rare, it 
is easy to lay out how it would occur. Consider a bead 
living at a Y-junction. We generically expect the bead 
to be moving, particularly since the growth of one of the 
connected strings is usually energetically favored. If this 
is the case, then it would only be a matter of time be- 
fore the bead at the junction would intersect any bead or 
anti-bead that lives on the strings that are "shrinking." 
When this collision occurs, the beads will interact, and 
we should expect their quantum numbers to add. 



the F-string is emitting massless axions. This gives it the 
long-range interaction typical of global strings. 

Global strings cannot break, so it is wrong to think of 
the F-string as "breaking" on a bead. What is happening, 
as we showed in Fig. (fT]) and @, is that two or more 
strings with different quantum numbers are joining on 
a bead. All strings still have long range interactions, so 
we can still detect their topological charge at infinity. 
It is therefore impossible for the F-strings to disappear 
by "breaking" on beads. Note that even though each 
string individually has a boundary (where the bead is), 
the whole cosmic string (the cosmic necklace with the 
beads) does not. 

Nevertheless something novel does happen at these 
junctions. The strings that join at beads have different 
charges, and charge must be conserved at the intersection 
of the strings. A D3-brane wrapping a 3-cycle in the ex- 
tra dimensions can do just that. Indeed, as we mentioned 
before, the KS throat supports M units of F3 = dCi flux 
over a 3-cycle in the extra dimensions, 
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FIG. 2: An interaction between an (M-l,0)-string and a (2,0)- 
string forming a (l,0)-string. Note that the Y-junctions now 
have beads on top of them. Because of that, it is no longer 
necessary for the tension to cancel exactly. If they do not the 
beads will accelerate. 



B. Stability and D3-brane as Beads 

The F-string is charged under the NS-NS 2-form B^. 
It has the standard Kalb-Ramond action 
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where m = l/2iva' is the F-string charge. The Kalb- 
Ramond action can be obtained from the action of a 
global Higgs string for distances much larger than the 
size of the core [24] • In 4D, the NS-NS 2-form is dual to 
an axion, so the Kalb-Ramond action just tells us that 



The action Q is normally invariant under gauge trans- 
formation of B2, 8B2 = d\\. However, since we now have 
a boundary (for each part of the string in the necklace), 
this is no longer true. The variation of the actions of all 
the strings ending on a given bead will add to 
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where the integration is over time. This must be can- 
celled by the D3-brane action. The relevant term in the 
action is the kinetic term for the gauge field F = dA 
living on the D3-brane, 
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where m = ( 2 ir) 3 a'' 2 * s ^ ne D3-brane charge. Varying this 
gives 
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By looking at the equation of motion for the gauge field, 
one finds that *(£?2 + 2na'F) = g s C2, leading to 



= -M3 / \\dC 2 , 
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This exactly cancels the contribution from the F-strings, 
showing that the beads are indeed wrapped D3-branes. 

As global strings, one might expect the F-strings (as 
well as the D-strings) to be unstable to domain wall 
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formation. The reason for this is simple: we have not 
yet seen a massless axion, and the presence of a mass- 
less scalar field is strongly constrained in cosmology (by 
e.g., big-bang nucleosynthesis tests). Therefore, the ax- 
ion must pick up a mass. We should expect to break the 
shift symmetry of the axion in doing so, creating a do- 
main wall in the process. The tension of the wall would 
then force the strings to shrink and disappear. 

One can avoid this disastrous situation if the axion is 
"swallowed" by a gauge field. This can be accomplished 
in string theory with a Green-Schwarz-like coupling be- 
tween B 2 and some gauge field [H, [25[ . The strings then 
become local instead of global, and the B2 fields decay ex- 
ponentially away from the strings, leaving no detectable 
long-range charge. Polchinksi has conjectured [26[ that, 
in string theory, all local strings of this type are unstable 
due to monopole formation. To see this in any particular 
case, we would have to look in detail at a complete string 
theory model where we could identify the monopole on 
which the string can break. These monopoles would be 
different from the beads we discuss in this paper. Unlike 
the domain wall instability, the monopole instability is a 
mild one: it is quite possible for the monopole mass to 
be so high that the strings are long-lived and metastable. 
For instance, the heterotic 50(32) string, which exhibits 
this kind of breaking, is effectively stable for all practical, 
cosmological purposes [26[ . 

Hence the question of stability is a model-dependent 
issue. For the remainder of this paper we will suppose 
that the strings are local and metastable, with a life- 
time longer than the age of the universe. Some progress 
has been made recently in identifying the monopoles that 
the strings would break on in the KS throat as a D5- 
brane configuration (27j . These monopoles have very 
large masses, in agreement with our working assump- 
tion of metastability. Furthermore, it shows that the 
beads we consider in this paper are not the same as the 
monopoles on which Polchinski's conjecture predicts the 
string would be unstable to breakage. These develop- 
ments suggest that there is much more work to be done 
in understanding the underlying topological structure of 
these strings. 



C. Binding 

In general, the tension formula for strings in a warped 
throat (Eq. [T]) leads to higher bindingenergies than in 
flat space. Ongoing numerical work [28| has suggested 
that vortex binding may occur efficiently in a gauge the- 
ory model. A recent analytical study of the kinematics of 
string binding [30] confirms that bound states can indeed 
grow relativistically. We thus retain the working assump- 
tion from Ref. 0] that strings bind instantaneously. The 
new structure in the (p, q) tension spectrum (Eq. []} leads 
to a different set of binding rules than is described in Ref. 
Q. The most important difference is that non-coprime 
(p, q) states are now bound (these are states for which p 



and q have a common factor; i.e., p — kN, q = IN). In 
that paper, the fact that these states were not bound was 
important in preventing UV-catastrophes in the network 
evolution (For explanation, see note at end of bibliogra- 
phy, [46| ). A related change is that (p, 0) states with 
p > 1 are also now permitted to exist. (0, q) states are 
not bound, however, and cannot form. Adding all these 
new states and removing an important string binding rule 
(the decay of non-coprime states) reopens the problem of 
cosmological UV catastrophes. We tested our model with 
these assumptions and found that the cosmological den- 
sity did continue to grow to late times. However, perhaps 
these new rules are too simple. Non-coprime string states 
still have much less binding energy than coprime states 
when p <C M. To account for this, we prevent some 
fraction of non-coprime states from forming by allowing 
them to "decay" immediately upon formation as before. 
Let us define 
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We use pdecay ag a branching ratio between the forma- 
tion of non-coprime states and their coprime daughter 
states. This immediately restores network scaling. 



III. NETWORK EVOLUTION 

To study the evolution of string networks in a curved 
background, we modified the network model described 
in Ref. [i[ to include the new features described above. 
Our review of this model will be terse; readers looking 
for more detail are referred to the original paper, where 
it is discussed at far greater length. The new features we 
have added to that model are as follows. 

• We have introduced a maximum value, M, for the 
quantum number p. Thus, we take (k + M, I) = 
(k,l), so that binding interactions that seem to 
form a p = k + M state actually populate the state 
with p = k. 

• When we attempt to form a p = k + M state, we 
must take into account bead formation. This intro- 
duces a tunneling suppression factor related to the 
bead mass; see Eqn. [7] 

• The set of allowed states now includes non-coprime 
(p, q) states as well as all (p, 0) states, since they 
are now bound. 

• Non-coprime states are allowed to decay, upon for- 
mation, with a probability pdecay^ 
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• The new tension formula, Eqn. [T] is used in the 
calculation of interaction probabilities. 

• We neglect bead dynamics in the evolution of the 
network. The justifications for this are detailed in 
the next section. 



A. The MTVOS model 

Following Ref. [i[ , we define an average string velocity, 
v, and single network string correlation length-scale, L. 
Each are shared by all strings in our network, regardless 
of their tension or number density. Using the cosmologi- 
cal scale factor, a, we can also introduce a conformal time 
by drj = vdt/a and a comoving string length, I = LI a. 
Denoting each string (p, q) state by a single Greek let- 
ter for compactness, we can further define a scaled string 
number density, 



2 a 2 p Q \/l - v 2 
l\ a = a n a = ■ 



and a string interaction probability P Q/ 3 7 (the probability 
that the interaction of a string of type (5 with a string 
of type 7 will produce a string of type a). Since it is a 
key point, we restate it for emphasis: in this model, the 
correlation length L is shared by all (p, q) states. It is not 
related to the physical distance between strings (The dis- 
tance between two strings of type a may be given roughly 
in terms of that state's number density; ~ N a 1 ) . We 
may then assemble a set of network evolution equations 
for evolving each of these string states, the overall corre- 
lation length, and the average string velocity. To connect 
conformal time to physical time, we assume a Friedmann 
evolution for a in a radiation or matter-dominated uni- 
verse. We take into account string loop formation (the 
terms proportional to C2 and P) ; string binding (propor- 
tional to F); string stretching (ci7j); Hubble friction for 
reducing v (ex. Ha); and acceleration of v due to string 
straightening (oc 
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with primes denoting derivatives with respect to confor- 
mal time; the constants c\ — 0.21 and c 2 — 0.18 are 
obtained by matching the scaling-era values of HL and 
v to results from old string network simulations, while 
the constants F and P are generally left as variable pa- 
rameters. The constant F indicates the overall strength 



of inter-string binding interactions; we take F = 1.0 for 
all our integrations. The constant P indicates the string 
intercommutation and loop production rate. Note that 
P is not directly related to the probability of string inter- 
commutation. It is better thought of as a loop formation 
efficiency parameter. Based on the results in Ref. [9|, 
we know that for non-binding networks (where F — > 0), 
we must take P = 0.28 to match to simulations. Hence, 
we simply assume this value for P in all of the present 
calculations. Varying P changes the overall number of 
strings present in the simulation - see 0] for discussion 
- but affects the results for beads minimally. 



B. New Binding Physics 

All of the new physics introduced by the new warped 
geometry under consideration are encoded in the string 
interaction probability term, P Q( 3 7 . First, we write down 
the general form of the formula for this, given in Ref 0] , 
taking /3 = [k, I) and 7 = (m, n) to be the two interacting 
string states: 
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where P + indicates a = (p,q) = (k + m,l + n) and P~ 
indicates a — (p,q) = (k — m,l — n). The critical angle 
that differentiates between the two type of binding is: 
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The use of this critical angle encodes the physics of the 
new tension formula. The next things to include are the 
suppression of string interactions that form beads and 
the periodicity of the string tension states. We include 
these effects in P in the following way: 



P 



( P +M,q)(k,i){m,n) — - r tunncling- r (p.<?)('M)('"' 



where ^tunneling = e 7 ' Mb ^ max is only included when 
the bound state being formed has quantum number p > 
M. 

Finally, we include the probability of decay for non- 
coprime states with small binding energy. Recall that 
this is necessary because in a warped geometry non- 
coprime states are slightly bound. The decay of non- 
coprime states was very important for the regulation of 
UV catastrophes in previous integrations of these equa- 
tions [9( . We originally thought that the periodicity of p 
might help to regulate these catastrophes; early numeri- 
cal experiments told us that this alone was not enough. 
To implement this decay, we take 



(Nk,Nl) 



N(k, I) with probability P deca y 
(Nk, Nl) with probability 1 - P deca y 
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In practice, we modify our network interaction term so 
that, when a = (Nk, Nl), 

P a __.fi- p dcca y ^ p i p dcca y p 

Note that the overall probability does not add up to one 
anymore. Effectively, the factors ^tunneling an< ^ P^ eca y 
reduce the overall factor of F for certain string interac- 
tions. Since these reductions only apply for some specific 
interactions, it is more convenient to include these ef- 
fects in P a /3-y than to make F also a function of the string 
quantum numbers. 

C. Network Results 

As we alluded to before, the immediate result of our 
early network integrations was that some sort of decay of 
the slightly-bound states is necessary for the networks to 
go to scaling. Without any decay mechanism, the number 
of allowed states appears to be too large, leading to a UV- 
catastrophe, albeit a mild one. The probabilistic decay of 
non-coprime states, however, very easily regulates these 
catastrophes. 

The chief results of this network model are quite similar 
to those of the simpler model studied in @. We should 
note that another group has proposed a related but more 
detailed model, which also returns very similar results 
(2_j |. In the numerical integrations presented here we 
allowed p to range between and M — 1 and q between 
and 3; we also always assume F = 1 and P = 0.28. 
Varying the maximum value of q had little effect on the 
results, while leaving out states with p < M can lead to 
serious errors. Our results: 

• Networks are still able to scale. 

• Networks reach scaling quite quickly for each value 
of M that we tested (see top panel of Fig. [3]). 
Larger values of M could not be calculated due to 
the large number of states involved. However, we 
expect to recover the results of Ref. [j| for networks 
with very large values of M. Since we would expect 
any novel network behavior to be most manifest at 
low values of M, it was sufficient to test just a few 
such values. The dependence on initial conditions, 
discussed below, is an instance of such a behavior. 

• For a given set of initial conditions, the overall den- 
sity in cosmic strings does not depend on M. For 
the simplest initial conditions (equal numbers of 
(0,1) and (1,0) strings being the only states ini- 
tially populated) we find, for g s = 1, 0. strings — 40 
for all runs (see Fig. [3]); for g s = 0.5, we find 
^strings — 70, where we have defined ^strings = 

^strings/ ((8/3)^6*^(0,1)) 

• String density is still very dependent on string ten- 
sion: we again find N a cx fi~ n , and n ~ 8 for 
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1.2 x 10 7 
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9.4 x 10 10 


20 


1.1 x 10 9 
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2.2 x 10 13 


5.1 x 10 12 


22 


2.1 x 10 10 


3.9 x 10 9 


1.2 x 10 15 


2.7 x 10 14 


15, p > M/2 


1.1 x 10 3 


35 


4.5 x 10 4 


100 


20, p > M/2 


624 


54 


1.4 x 10 3 


170 



TABLE I: The relative populations of the lowest-tension p = 1 
and p = M — 1 states. The first four rows were calculated 
under the standard initial conditions, for which only p < M/2 
states are initially populated. The last two are calculated 
under initial conditions where all states, including p > M/2, 
are initially populated. Note that = fJ*(M-l,q)' In au 

cases, the (M — 1, ±1) state had a larger population than the 
(M - 1,0) state. 

standard initial conditions. However, this relation 
is now more an envelope function than a predic- 
tive function: the sinusoidal shape of the p tension 
leads to greatly-enhanced structure in the string 
spectrum (see Fig. 2]). For instance, for M=15, the 
tension of the (1,0) strings and the (14,0) strings 
are identical; but, for our standard initial condi- 
tions, -V( li0 )/-V( 14)0 ) = 2.5 x 10 6 . See table [Q for 
more such comparisons. 

• The complexity of the string tension spectrum leads 
to a dependance on initial conditions in £l s trmgs- 
Essentially, the sinusoidal shape of the tension 
formula for p strings suppresses the formation of 
p > M/2 states if the initial distribution of string 
states is confined to p < M/2. However, if the 
p > M/2 states have similar initial populations to 
the p < M/2 states, then the scaling value of fl is 
approximately doubled. For instance, for g 8 = 1, 
the scaling value of jumps to 95 when states with 
p > M/2 are initialized. See Fig. [5] for an example 
of this effect. 

• The tension spectrum is also bimodal in initial- 
condition response. The results pictured in Fig. 
2] show the tension spectrum for initial conditions 
where the p > M/2 states are initially populated. 
In that case, the envelope-function for number den- 
sity versus tension is much less steep: n ~ 3 or 4. 



D. Other Considerations 

We also considered including an interaction term that 
took into account the strings' kinetic energy in deter- 
mining whether a particular binding interaction would 
occur. We wrote down a suppression term similar to 
our expression for pdecay^ exce pt with a denominator 
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FIG. 3: The bottom panel shows the evolution of D. cs = 
f2os/((8/3)7rG/i(o,i)) for various values of M, taking the string 
coupling g s = 0.5. The top panel shows the rate of change 
in the comoving number density Nr/ 1 ; in the scaling regime, 
dlog Nrj 1 /dlogrj = 0. 
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FIG. 4: The final scaling-era spectra for M = 20 taking the 
string coupling g s == 1.0, with -/V( ,i) normalized to unity and 
the other number densities scaled accordingly. The black stars 
are for the standard initial conditions; the red triangles are 
for initial conditions in which the p > M/2 states are initially 
populated. String spectra for other values of M are similar. 
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FIG. 5: The evolution of Q cs — f2 cs /((8/3)7rG/i(o,i)) for vari- 
ous initial conditions, taking the string coupling g s = 1.0 and 
M = 20. The black, solid line: only (0,1) and (1,0) states are 
initially populated. The blue, long-dashed line: all states with 
p < M/4 initially populated. The red, dashed line: only (0,1) 
and (1,0) states, but with the initial populations quadrupled 
in size. The black, dash-dotted line: all states with p < M/2 
populated. The blue, dot-dashed line: only the (0,1), (1,0), 
and (M-1,0) states are initially populated. The magenta long- 
dashed/dotted line: every allowed (p,q) state is begun with 
a high initial population. We see that there are two scaling 
values of Q: ~ 40, for cases where the initial conditions have 
only the low values of p < M/2 initially populated; and ~ 95 
when the states with p > M/2 are initially populated. 



based on the interacting strings' kinetic energy. Evolv- 
ing networks with this suppression term did regulate the 
most obvious uncontrolled growth of cosmic string net- 
work density, even when we did not allow the decay of 
non-coprime states. However, the networks did not enter 
a scaling regime: instead, they continued to evolve to late 
times, not becoming dominant but with an always slowly 
increasing cosmological density. The reason for this is 
that kinetically suppressing string interactions slightly 
favors the growth of high-tension number density. Heavy 
strings have a great deal of kinetic energy. If strings 
with high kinetic energy cannot easily undergo binding 
interaction, heavy strings become isolated; their effective 
interaction probability decreases, F — > 0. Since phase 
space and binding energy effects both tend to decrease 
high-tension-state number density, we get a competition. 
For values of M < 15, this competition leads to late-time 
network evolution, even when the decay of non-coprime 
states is allowed to proceed. The MTVOS model, as for- 
mulated, is difficult to evolve to late times, since steps 
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in conformal time drj oc a -1 . It may be worth reformu- 
lating the equations to study this late-time network de- 
velopment, which may provide additional string network 
phenomenology that would be interesting to study fur- 
ther. For the time being, though, we have simply omitted 
any kinetic energy-based suppression from our model. 

While this paper was being prepared, Copeland, Kib- 
ble, and Steer [30j] published a detailed analysis of kine- 
matical constraints on string binding, building on some 
results which they had previously published [3l|. Their 
formulae for kinematical determination of whether bind- 
ing proceeds are derived in the Nambu-Goto string limit 
for several classes of string binding configurations. Since 
our network model needs binding probability formulae 
that cover all possible binding events, we were not able 
to apply their formulae to our model. However, two of 
their chief results - that lighter strings are preferred by 
kinematics and that string binding, when it occurs, pro- 
ceeds relativistically - are in line with the assumptions 
and conclusions of our work. Any future development of 
the MTVOS network equations should take into account 
the constraints that Copeland et al. have derived, in- 
cluding their finding that binding networks may exhibit 
smaller r.m.s. string velocities. Ultimately, a fusion of 
their approach with the approach taken here could pro- 
vide a robust framework for the quasi-analytical study of 
binding string networks. 



E. Populating Low- Tension, High-p States 

In Ref. [9[, it was proposed that the main observable 
difference between a network of (p,q)-strings and a sin- 
gle string network is a factor of 3 in fl cs . The reasoning 
behind this is that the network reaches scaling with sig- 
nificant populations in only its lowest-lying tension states 
- those not liable to decay. For the (p, g)-strings network 
in flat space, these are the (1, 0), (0, 1) and (1, 1) states - 
a total of 3 states leading to a factor of 3 in energy den- 
sity as compared to a single-species string network. For 
a general network with "m" types of strings, the density 
of strings should, then, be 2 m — 1 times the results for a 
single string network. 

In our case we have two additional low tension states, 
(M — 1,0) and (M — 1,1). This, too, leads to an enhanced 
density of strings. However, this enhancement is appre- 
ciable only if these states are populated initially, rather 
than being formed only through dynamics. When these 
sorts of initial conditions are considered, we can see from 
Fig. [5] that the overall network density is approximately 
doubled, just as this rule of thumb predicts. 

Initially populating high quantum number states is, 
however, probably unrealistic. Also, these initial condi- 
tions may already be excluded, since highly populated 
(M — 1, n) states can lead to problems with overproduc- 
tion of beads. We will discuss this concern more below. 



IV. BEADS ON A STRING 

As we have argued before, whenever two strings meet 
such that the absolute value of the sum of their p quan- 
tum numbers exceeds M, a bead must be formed for the 
interaction to proceed. The beads that are formed do 
not exist on their own. They are constrained to live on 
their confined flux tubes - to be threaded by the string 
network on which they were formed. Since the 0-string 
and the M-string have zero tension, the beads are never 
a monopole-like "end" of a string. Instead, they must 
function as "beads" on the string, always connecting two 
different kinds of string 13]. As mentioned previously, 
bead production is a tunneling process, and therefore it 
can be totally suppressed if the bead mass is very large. 
When M is not too large, though bead production can- 
not be neglected. We might worry that beads, as point- 
like sources of energy-density, could cause serious prob- 
lems for the cosmology of this model. Indeed we know 
that ordinary monopoles can overdose the universe, since 
their densities decrease only like 1/a 3 and hence grow 
rapidly with respect to the radiation component of the 
universe until they are dominant. In our case, since these 
monopoles - or beads - are stuck on strings, we expect 
a very different behavior. A useful parameter to study 
this is r — where d is the average distance between 
beads 

For this first analysis we shall consider only a single 
cosmic string tension. The average mass per unit length 
of such a necklace is (r + 1)^. For r <C 1, the beads 
contribute a negligible amount to the mass density on a 
string and we can suppose that the beads are simply car- 
ried around and remain cosmologically innocuous. For 
r > 1, the beads are the dominant contribution to the 
mass density and the network should look like a gas of 
monopoles. This is the cosmologically disastrous case. 
Intermediate values of r (r ~ 0(1)) are interesting, be- 
cause they may not be ruled out, yet would imply a higher 
effective string tension and energy density for the string 
network. 

In the cosmological model presented in this paper we 
start with a network of (jp, q) strings. These strings pro- 
duce beads only through interactions - there are no pri- 
mordial beads. Hence we necessarily start with rC 1. 
What we need to calculate is how this dimensionless num- 
ber evolves with the development of the network. A sim- 
ple model due to Berezinsky and Vilenkin describes the 
salient features 1231. 



A. Evolution in the Scaling Regime 

To begin simply, we will assume that a cosmic super- 
necklace network has reached a scaling solution without 
having formed too many beads, so that it enters the scal- 
ing era with r -C 1. We can then study the time evolution 
of r by looking at energy conservation for a particular 
(closed) necklace of length L = Nd — N * , where N is 
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the total number of beads on the necklace, which we take 
to be a constant. The size of this loop- necklace can be 
greater than the Hubble length; if so, the closed necklace 
will appear to be an infinitely long string in any particu- 
lar Hubble volume. In the spirit of Ref. [22j], we can then 
write an energy balance equation for this closed necklace: 



Gravitational Radiation 



E 



-PV - E, 



(20) 



where E is the total energy of the necklace, P is the 
effective pressure and E g is the energy loss by gravita- 
tional radiation from small scale wiggles. This approach 
neglects several things. 

• The production and annihilation of beads. We will 
return to this point later and discuss the effect of in- 
cluding bead production. Annihilation could be an 
important effect when r approaches 1 while unim- 
portant in the scaling regime which we discuss be- 
low, where the inter-bead distance is of order the 
Hubble scale. However, it might be important in 
the transient regime. Given the uncertainties, we 
did not include annihilation in our analysis. Inclu- 
sion of annihilation will only aid in preventing bead 
domination of the network's dynamics. 

• Energy lost to gravitational radiation from moving 
beads. Assuming r -C 1 and that no average force 
acts on the beads should keep this contribution neg- 
ligible. 

• Loop formation. We will come back and include 
the effect of loop formation later. Loop formation 
can be very important. In fact, we will argue that 
it leads to a scaling value for r; but to get to that 
point we need to understand the evolution of r in 
the absence of loops first. 

• The physics of Y-j unctions, which can "capture" 
beads and transport them around as we have dis- 
cussed briefly and will touch on again below. We 
assume this will not change the physics too much, 
but we have not yet been able to quantify this as- 
sumption. 

• Finally, the fact that the work term — PV is really 
only thermodynamically well-defined for an ensem- 
ble of cosmic necklaces. Nevertheless, we can con- 
nect this averaged quantity to a particular loop by 
considering the work done per cosmic necklace in 
the scaling regime. This is, of course, a rough ap- 
proximation. 



Eq. (|20[) contains two terms through which the energy of 
a necklace can change. The first term is the work done 
by the Hubble expansion. Depending on the sign of the 
pressure, this can cause either a gain or a loss of necklace 
energy. The second term is the loss of energy due to 
gravitational radiation. Below, we will assume that most 
of the gravitational radiation is coming from the small 
wiggles on the string. 



The power emitted by the wiggles on a string is ap- 
proximately given by [24j | 



dP 



H) 2 w 



(21) 



where Co is the average frequency of the wiggle and jl is the 
effective tension of the string. Estimates of the effective 
tension of a wiggly string vary, from 2/i in standard lore 
[24j to 1.1/x in more recent analytical calculations [32| . 
The precise value turns out to have very little effect on 
our result. Nevertheless, we leave it as a free parameter 
and write jj = (3fi. We take the average wavelength of 
the wiggles to set the size of the loops Zi oop = at and 
therefore Co — 2% /(at). With all this we find 
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dl 
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(22) 



where L — Nd and k. 



_ 4tt 2 Gm(/3-1) 



2. Hubble Stretching 

To calculate the Hubble stretching of a single cosmic 
necklace, we first calculate the Hubble work on the net- 
work as a whole during scaling. Using this result, we can 
approximate 



(--PVOl-aecklace « (~PV) 



ncckla 



total " 



E, 



total 



(23) 



Assuming that r <C 1, we can suppose that the pressure 
is mostly due to a gas of cosmic strings. We then have 



P 



1) 



(24) 



where p s <~ ( r + 2 1 )^ t during scaling and £ is the correlation 
length in the scaling regime. For a network of cosmic 
necklaces, we expect £ ~ (r + l)it Taking the to- 

tal volume to be a comoving volume characterized by a 
length scale, L, we get 



(-PV) total 



^Vk s 



(25) 



where V/t; 2 is the total length of strings in the comoving 
volume V, k s = (1 — 2v 2 )v and we took a ~ t v . The 
total energy of comic necklaces in a comoving volume V 
is E = (r + l)/iV/£ 2 and the energy of the specific cosmic 
necklace we are looking at is E = (r + l)/j,L. Hence the 
Hubble work for 1-necklace of length L is 



(-PV)x- 



NM bKs NM b K s 



-necklace 



rt 



rt 



(26) 
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This is of course an approximation, but it has the correct 
asymptotic behavior, i.e., it grows for smaller value of r. 

Plugging everything in the energy balance equation 
(Eqn. [20]), we obtain 



t 



f 



(27) 



Using the scaling value for v 2 from simulations [H, [H| , 
we find that n s is 0.07 and 0.14, in the radiation and 
matter eras, respectively. This gives a negative cosmic 
string pressure, so the strings are stretched by the Hub- 
ble expansion. Thus, the energy balance equation for 
the necklaces expresses a competition between Hubble 
stretching and shortening due to gravitational radiation. 
Note that at r = and/or t — 0, the preceding equa- 
tion fails. This is not a serious defect, just a limitation 
caused by our way of formulating the problem. This for- 
mulation is convenient for us, and since we will always 
consider r > and can take our initial time (time at 
the end of inflation) to be non-zero, it will not cause any 
problems for our analysis. 



3. Solution and Interpretation 

The solution to this equation is r(t) ~ t Ks ~ K '. This 
result can be understood as follows. The stretching of 
the string (k s ) tends to decrease r, whereas the radiation 
from the small scale structure (k 9 ) tends to shorten the 
string and increase r. 

In Ref. [13], they assumed that K g ~ 1, which follows 
from the assumption that a ~ G/i. This predicts an r 
that grows until it becomes a value of order 1, at which 
time the various approximations done in this calculation 
break down. If, on the other hand, we take at face value 
the recent results from [35ll3r| that give a ~ 0.1, we then 
find k 9 < 1 and r decreases with time. We should also 
note that other recent expanding space simulations sim- 
ilarly find considerable loop production at length-scales 
much larger than the gravitational back-reaction scale 
[37L [38j . In this scenario, as a grows, the small-scale 
structure gets bigger, and we thus get less gravitational 
radiation for a given length of string. The stretching then 
wins over the shortening and we find that the beads will 
get farther and farther apart. 

This conclusion is also consistent with the recent an- 
alytical results of Ref. 32], who found that stretching 
is more important than what was previously anticipated 
from simulations. Given these successes, we claim that 
this approach, though clearly very simplified, captures 
the essence of the problem. On a given string we should 
expect the beads to get farther and farther apart if the 
string stretching from Hubble expansion is more impor- 
tant than the string shortening from gravitational radia- 
tion. Of course, the value of a is still a subject of debate. 
It is hard completely to determine the status of the prob- 
lem without a more in-depth analysis. A reader wishing 



to pursue these questions more will find a good starting 
point for such a detailed analysis in Ref. [391 ] - 

We can also take an alternative approach and try to 
determine a range of a values such that the stretching 
dominates over the shortening. From the numbers given 
above, this range is given by 



a> 10 3 Gn{(3-1) 2 . 

With G\i ~ 10~ 7 this gives roughly 

a > 10~ 4 for p = 2, or 
a > 10~ 6 for /3 = 1.1 



(28) 



(29) 



These values for a lie in the "middle ground" between 
the various claims in the literature, so finding the real 
value of a is critical for determining the fate of beads on 
cosmic necklaces. 



4- Including Production 

One can include the dynamical production of beads in 
the previous analysis by allowing N to be a function of 
time. Since we want to find conditions for which beads 
are subdominant, we take r <C 1 as an ansatz. Using 
this, we find: 



t 



t 



N 



(30) 



When we reach the scaling regime, a certain number 
of beads have been formed. The network has settled into 
a stable state. The great majority of the strings are in 
the lowest energy states: (1,0), (0,1) and (1,1) (and, 
possibly, (M —1,0) and (M — 1,1), depending on initial 
conditions). In the former case, it is only rare interac- 
tions that produce beads. We can estimate this rate in 
the scaling regime very simply. The number of collision 
per unit of time is ~ 1/t (for strings going at speed v ~ 1 
and separated by L ~ t) , and the number of beads formed 
per unit of time is the number of such collisions that form 
beads. Looking at Table IT] we can guess that the domi- 
nant bead- forming interactions will be between (1, n) and 
(M — l,n) strings, where n = Q, 1. The smaller popula- 
tion of the (M — 1,0) states will be the first source of 
bead production suppression. As a proxy for this effect, 
let us define, 



R : 



-W(M-1,0) + N (M-1,1) 
N (l,0) 



(note that R is related to the inverse of the ratios as 
given in Table [TJ) . The second source of suppression is 
the tunneling suppression. Together, these imply a rate 
of production given by something like 



N 



R 
— ( 

t 

R 



-■Ptu= 



-M 



(31) 
(32) 
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There is no clear relationship between M and R, but for 
most scenarios we should expect R to decrease for in- 
creasing M; it will certainly never increase drastically. 
Thus, tunneling suppression forces this rate to zero for 
large values of M. The rate also decreases at late times, 
when strings gets farther and farther away from each 
other and collisions become rarer. 

Still, N(t) ~ lnt, so the number of beads always grows 
- just very slowly. Thus, we can make the term N/N 
negligible in (|30| because of its very small coefficient. 

Although production is negligible in the scaling regime, 
it will be important during the transient regime. One 
might reasonably fear that the presence of too many large 
quantum number strings at the onset of the network evo- 
lution might create too many beads. 



B. Loops and the Scaling of r 

In this section we discuss the effect of loop formation 
in the scaling regime. Loops can have a very important 
effect on the value of r. In fact, loop formation should 
promote a scaling solution in some cases. To see this, we 
first observe that each loop that is formed removes an 
amount I of string length. When a loop is formed, there 
may or may not be beads that are trapped on the loop. If 
there are no beads on the loop, then the loop formation 
will decrease the distance between beads, increasing r. 
If there is a bead and a anti-bead on the loop-forming 
portion of the string and the average interbead distance 
is half the average size of loops, then r will be unaffected. 

If r would decrease from stretching in the absence of 
loops, adding loop formation can either slow down the 
decrease of r or cause it to cease decreasing altogether. 
The latter is the interesting case. 

In necklaces with sparse beads, loop formation removes 
length between beads. As this process goes on, the av- 
erage spacing between beads will eventually become of 
order half the average loop size, I — at. At this time, 
we will expect nearly every loop to contain a bead and 
an anti-bead. This should slow or stop the evolution of 
r from loop formation. This argument suggests a scal- 
ing value for r of 2Mb/ (fj£), though the contingencies of 
a fully realized network will doubtlessly complicate this 
simple picture. Nonetheless, this mechanism seems to us 
robust, and should enforce a lower limit on r for cosmic 
necklaces that can form loops. 

Note that if r grows in the scaling regime, then loop 
formation will only enhance that effect and there will be 
no scaling in r during an epoch of rapid bead production. 
Hence we find that if a < 10~ 4 (from section PV A 3|) . 
then r will grow. If this is the case, then cosmological 
considerations will force M to be huge so as to suppress 
bead formation (see next subsection). 

On the other hand, if a > 10~ 4 and Hubble stretching 
is tending to spread beads apart, then we find that r 
scales to the following range of values (taking a ~ 0.1 as 



an upper bound): 



10M 3/ 



2 y/2^b H 
37r 3 / 2 m s 



< r < 



10 4 M 3/2 



V2{hb H 
37r 3 / 2 m s 



(33) 



Interestingly, the value of r grows for large values of M. 
Beads are potentially observable (though our analysis be- 
comes less trustworthy) for values of order 1. This would 
require very large values for M, though, since H/m s is 
such a small number after inflation. Unfortunately, since 
H decreases with time, it seems unlikely that the beads 
will ever be observable. A value for r of order 1 at the last 
scattering surface (with potential observable effects in the 
CMB) would mean a value of r bigger than 1 at earlier 
times, where our analysis can no longer be trusted and 
that very likely would lead to other cosmological prob- 
lems. Alternatively, for the network to have r of order 
1 at the earliest time where we could hope for a lensing 
event to be observed - at, say, a redshift of 1 - would 
mean a value of r of about 10 4 at the surface of last scat- 
tering, leading to disastrous cosmological consequences. 
We also note that, dynamically, large values of M ac- 
tually suppress the creation of beads: it might thus be 
impossible for any dynamically-evolving network to form 
enough beads to realize this scaling solution for r. 



C. Evolution in the Transient Regime 

It is much more difficult to analyze cosmic necklaces 
in the transient regime. In our network evolution, we 
have assumed that bead dynamics can be neglected. Al- 
though this is definitely true in certain cases (when M 
is very large, for example), we might expect that certain 
initial conditions, together with certain values of the pa- 
rameters, might lead to a large production of beads. 

One thing we can say for certain: in the transient 
regime, r initially grows. On the one hand, a is expected 
to be of order G[i in the transient regime [33], Hj, l35j ] . 
leading to more gravitational radiation production per 
unit length. This should cause the beads to get closer 
to each other. Furthermore, we expect that the Hubble 
stretching is going to be less efficient before we reach the 
scaling regime (since the string-gas pressure goes from 
positive to negative during the transient regime). Fi- 
nally, a large initial density of strings or a significant 
initial population of high quantum number states can 
lead to considerable bead production in the initial stages 
of the evolution, when higher physical density will lead 
to far more string interactions and thus bead formation 
opportunities. 

To get a feeling for how critical this situation is, we 
can do some simple approximate calculations. To begin, 
we neglect bead production and loop formation. We will 
also assume that Hubble stretching is negligible. Taking 
k s = 0, we find r ~ t Kg = r Q t i ' nGtJ '^ a . Since we know 
that the scale factor "a" grows by at least 3 orders of 
magnitude before we reach scaling (which corresponds to 
6 orders of magnitude in time during the radiation era), 
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we can estimate how small tq needs to be assuming r is 
of order 1 when we reach scaling. We find 

r ~ KT 6 ^)^/" ~ 1(T 107 , (34) 

where we have assumed the most extreme case: G\x ~ 
10~ 7 and a ~ 10~ 12 in the transient regime. More re- 
alistic values for a will result in much larger and more 
reasonable upper limits for r$. From this calculation, 
we see that we need to start with essentially zero beads. 
Supposing that large quantum number (p, q) states are 
initially populated, then the only thing available for sup- 
pressing bead production is the tunneling coefficient. We 
get that r ~ N(t) ~ e~ M ln(£) from production. This 
means that ro grows like a log from production alone. 
Thus, we would need e~ M to be of order 10 -10 , which 
implies M ^> 10 3 5 . Again, this is the worst case sce- 
nario. A better situation occurs if a ~ G/i. In this case, 
we find M ^> 10 2 . Let us emphasize that this analysis 
assumes that the high quantum states are present ini- 
tially and that the only suppression is from tunneling. 
On the other hand, we are still neglecting loop formation 
as well as intercommutation, each of which increase r and 
which will, thus, force us to require an even larger value 
for M. Hence, if we have initial conditions in which high 
quantum number states are populated, we expect that 
M needs to be even larger than our previous estimates 
for beads not to dominate. It is hard to accurately esti- 
mate this quantity, since we do not know how effective 
loop formation will be. We conjecture that M of order 
of a thousand should be sufficient to keep the beads from 
dominating in the transient regime. In the less dangerous 
case - when high quantum number states are not initially 
populated - these requirements would be somewhat re- 
laxed, though a much more detailed analysis would be 
necessary to quantify by how much. 

Finally, we expect that the number of beads produced 
depends on the initial density and spectrum of the string 
network the same way the final scaling density did. This 
is very interesting, as it puts an indirect bound on quanti- 
ties - the initial string populations - that are usually left 
totally unconstrained in standard string network evolu- 
tion. Indeed, cosmic string networks usually reach iden- 
tical scaling solutions for any initial density of cosmic 
strings. On the contrary, we expect this not to be gener- 
ically true for cosmic necklaces. We do not believe this 
bi-modal scaling structure is an artifact of a poor state- 
ment of the problem, but in fact represents two physically 
distinct, initial-condition-selected scaling solutions. 

The Kibble mechanism predicts roughly one cosmic 
string per Hubble area at the time of formation. For 
cosmic superstrings, there are some indications that this 
number may be much higher, of order of the string scale 
[40| . Though such networks can reach scaling, the very 
large number of interactions they undergo at early times 
could lead to dangerously enhanced bead creation. Fur- 
thermore, it remains to be seen whether the same kind 
of analysis will suggest that the high-quantum-number 



states will be initially populated. Clearly this is yet an- 
other issue that needs to be explored further. 

In summary, we expect that the transient regime puts 
important constraints on both the value of M and the 
initial distribution of strings. However, we cannot fully 
quantify these bounds without a more detailed analysis. 

D. Multi-tension String Networks with Beads 

In addition to small scale structure versus stretching, 
a binding, multi-tension string network brings in other 
possible bead interactions. If an unzipping Y junction 
meets a bead, for instance, it will "grab" it. There will 
subsequently be three strings attached to the bead. Since 
the outgoing or ingoing charge will remain the same, we 
can always tell if it is a bead or an anti-bead. If two 
Y-junctions meet and the net bead number is non-zero, 
then one could create an X-junction. This process could 
in principle lead to junctions with any number of legs 
(perhaps with strange gravitational lensing properties cf. 
|4l|). However, for reasons similar to those given above, 
we expect the formation of such bead-bound junctions 
to be rare. Furthermore any junctions with more than 
three strings attached to it should be unstable to decay 
into multiple Y-junctions. 



V. CONCLUSION 

Multi-tension cosmic string networks that are formed 
and live in phcnomcnologically viable warped throats ex- 
hibit cosmologically safe scaling behavior under the fol- 
lowing assumptions: (1) string binding proceeds as in 
Ref. Q (2) Non-coprime states sometimes dynamically 
decay, even though they have a nonzero binding energy. 
We find n strmgs ~ 40(70) {8n/3)Gfi F for g s = 1.0(0.5) in 
our model, assuming the simplest initial conditions. The 
warped geometry causes these strings to have a tension 
that is periodic in i^-string charge, p. We find that this 
property causes the results of network evolution to be 
dependent on initial conditions: if strings with p < M/2 
are the only states initially populated, we find the above- 
stated results; but if we start with even a small initial 
population of p > M/2 strings, the scaling value of is 
more than doubled. 

Additionally, the formation of any bound states that 
"pass through" the maximum value for p leads to the for- 
mation of beads on the string. This causes the string net- 
work to become decorated with D-brane beads: a cosmic 
supernecklace. These beads, or monopoles, are cosmolog- 
ically safe if the average loop size is larger than 10 3 G/ii 
in the scaling regime, so long as their production in the 
transient regime is sufficiently suppressed. Including an- 
nihilation in the analysis could potentially relax this con- 
straint even further. In the transient regime, we have fur- 
ther assumed that the number of beads is small enough 
that it does not affect the network evolution. This con- 
dition is dependent on the mass of the beads (related 



14 



to the flux number M) as well as the initial distribu- 
tion of string populations. We expect these cosmological 
constraints to give direct bounds on M, though more de- 
tailed analysis will be needed to quantify those bounds. 
Finally, we argued that the inclusion of loop formation 
could lead to a scaling solution of the interbead distance. 
It might therefore be possible for beads to have observ- 
able effects, though realizing such a scenario appears to 
require fine-tuning. In any event, the beads' potential for 
creating distinctive observational phenomena should be 
more thoroughly explored. 

Looking forward, the methods discussed in this paper 
should allow us to study even more general multi-tension 
string networks. It was conjectured in Ref. [ll| that the 
general (p,q) string tension in a flux background should 
be given by 

T M ^mWg(q) 2 + f(p) 2 - (35) 

With appropriate parametrization, it should be possible 
to study the behavior of networks with general g(q) and 
f(q) using the techniques outlined here. As cosmologi- 
cal tests of cosmic strings become more stringent, such 
schemes may allow us to place useful limits on the struc- 
ture of the string tension spectrum. This will give string 
theory model builders guidance in constructing cosmolog- 
ically viable, phenomenologically distinctive theories for 
cosmic string formation. While the "square-root" tension 
spectrum may not be a distinctive signature of string the- 
ory, it is rather hard to reproduce it in standard effective 
field theory framework (see Ref. [42| for an attempt). 
Saffin has proposed an effective field theory model with 



two U(l) fields to mimic (p,q) bindi ng [43| , an approach 
anticipated by Carroll and Trodden [44| . Extending such 
models to produce periodic tension formulae is the next 
challenge. 

Cosmic superstrings remain one of the most compelling 
potential observational windows into physics at the string 
scale. They appear to be cosmologically safe; they can be 
created and maintained as metastable in realistic string 
theory models; and, as we have shown, their phenomenol- 
ogy is rich and varied. If they do exist and can be ob- 
served, they might provide unparalleled observational ac- 
cess to the physics of string theory. 
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[45] As stated, the tension formula for the (p,q)-string does 
not include the contribution from the fluxes it sources. 
This is usually a small, second order contribution, but 
in the case of the (M,q)-string - where the first order ef- 
fect vanishes - it is the dominant one. This contribution 
slightly breaks the periodicity of the string tension. This 
means that the M-string is actually a physical object and 
therefore while the process in Fig. 1 is still possible, one 
could also possibly form a junction with the very low ten- 
sion (T ~ (Gfi) 2 ) M-string without forming beads. Given 
the tiny phase space for junctions yielding extremely low 
tension string states, we did not include this type of in- 
teractions in our analysis. However, since this interaction 
will only further limit the production of beads, and thus 
make them even safer cosmologically, we do not consider 
this to be an important omission. We also note that such 
interactions should not appreciably change the network 
dynamics. Nonetheless, this deserves further study. 

[46] By UV-catastrophes, we mean the fact that there is no 
upper bound to string tension. Thus, we can imagine 
cases where strings continually bind together, increasing 
the quantum number with each interaction, causing the 
tension to evolve to ever-higher values. This sort of ever- 
upward evolution would be cosmologically disastrous, as 
very-high-tension strings would come to dominate the 
universe's energy density. 



